Analytic properties of unitarization schemes 

O. V. Selyugin *) 
BLTPh, Joint Institute for Nuclear Research, Dubna, Russia 

J.-R. CUDELL*) 

Institut de Physique, Bat. B5a, Universite de Liege, Sart Tilman, B4000 Liege, Belgium 
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impact parameter representation at high energies. Different unitarization procedures and 
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1 Introduction 

Saturation is now a very popular term, but it has a very wide meaning. It 
includes shadowing and antishadowing processes, the saturation of the Froissart- 
Martin bound [I], gluon merging at small x and geometrical scaling, unitarization 
effects, etc. In the S'-matrix language, saturation means that we reach the maximum 
possible scattering. This happens for a specific value of the impact parameter, and 
is directly connected with the unitary property of the scattering amplitude. It is 
that last meaning that we shall use in our work. 

The most important results on the energy dependence of diffractive hadronic 
scattering were obtained from first principles (analyticity, unitarity and Lorentz 
invariance), which lead to specific analytic forms for the scattering amplitude as 
a function of its kinematical parameters — s, t, and u. Analytic S'-matrix theory 
relates the high-energy behaviour of hadronic scattering to the singularities of the 
scattering amplitude in the complex angular momentum plane. One of the impor- 
tant theorems is the Froissart-Martin bound, which states that the high-energy 
cross section for the scattering of hadrons is limited by 

aZ? x = ^log 2 (-), (1) 

where s$ is a scale factor and fi the lightest hadron mass (i.e. the pion mass). As 
the coefficient in front of the logarithm is very large, "saturation of the Froissart- 
Martin bound" usually refers to an energy dependence of the total cross section 
rising as log 2 s rather than to the actual limit for the total cross section. 

Experimental data reveal that total cross sections grow with energy. This means 
that the leading contribution in the high-energy limit is given by the rightmost 
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singularity in the complex- j plane, the pomeron, with intercept exceeding unity. In 
the framework of perturbative QCD, the intercept is expected to exceed unity by an 
amount proportional to a s [2J. At leading- log s, one obtains a rightmost singularity 
at J — 1 = 12 log 2(a s /ir). Such a singularity — the hard pomeron — seems to be 
present in inelastic diffractive processes [3] and may be present in soft scattering, 
where it appears as a simple pole with intercept « 1.4 for energies smaller 
than 100 GeV [H[S]. In this case, the Froissart-Martin bound is soon violated, and 
the BDL-regime may appear at relatively low energies. The effect of saturation on 
the growth of the total cross section is however far from clear, because it involves 
processes at very small x and non-perturbative effects. 

Saturation in the framework of perturbative QCD may be connected with the 
growth of the gluon density at small x, which must be bounded by non-linear effects. 
In principle, the gluon density may have a maximum which will be reached before 
the BDL (see e.g. [5]). In fact, the non-linear corrections added to the BFKL equa- 
tion lead to a saturation regime at relatively large impact parameters, and are not 
directly related with the BDL, which first occurs at small impact parameters. One 
must note that, even in perturbative QCD, the description of saturation at large 
impact parameters is problematic because of non-perturbative effects, connected 
with confinement [7j. 



2 Unitarity: schemes and equations 



Unitarity of the scattering matrix 

SS+ = 1. (2) 

is most suitably studied in the impact parameter representation [5] as it is equivalent 
(at high energy) to a decomposition in partial-wave amplitudes. 

We thus take the scattering amplitude in the impact parameter representation 

/>oo 

T(s,t)= / bdbJ (bA)G(b,s). (3) 
Jo 

with 

ImG(b,s) < 1. (4) 

and 

ImG(s, b) = [Im(G(s, b)] 2 + [ReG(s, b)} 2 . (5) 

As energy grows, the scattering amplitude in the impact parameter representa- 
tion can exceed the unitarity bound for a range of impact parameters b < 6j. To 
restore unitarity, there are several different schemes. Two of them are based on the 
solution of the unitarity equation [9]. The first and most popular scheme seems to 
be the eikonal representation. 



T(s,t) = i / bdbJ (bq)(l - exp(-xM)). (6) 
Jo 
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where the eikonal function x(s,b), which in the non-relativistic regime is equal to 
the scattering potential, is often taken as the Born scattering amplitude. 

dN/dy 




0.2 0.4 0.6 



Fig. 1. The derivative of the non-linear equations corresponding to the eikonal (hard 
line), the rescaled [/-matrix (dashed line) and the hyperbolic tangent (long dashed line) 

unitarization schemes. 

The second unitarization scheme which we shall consider is the [/-matrix scheme 
[101 lllj , which we rescale so that the unitarized amplitude becomes asymptotically 
i as s — > oo: 

T( S ,i) = r bdbJ (bq) ( X{S f h j ) ■ (7) 
Jo V 1 -iX{s,b) J 

with t = —q 2 - If we assume that the Born scattering amplitude factorises as 
f(b)s 1+A , we can rewrite this scheme as the solution of the non-linear equation 

^L = AN[1-N}). (8) 

We find that eikonal representation also corresponds the non-linear equation like 
above [T^] . 

dN 

^- = -Alog(l-AT)(i-AO. (9) 
dy 

We can also obtain the above schemes by assuming that dN/dNsom = D(N), in 
which case 

dN = dNBorn (1Q) 
dy dy 

We can also easily obtain new unitarization schemes [13j . such as 

T(s,t)= bdbJ (bq) tanh[x(s, &)])• (11) 
Jo 

which corresponds to the non-linear equation in the form 

dy dy 
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At first sight, these three schemes are quite different. However we can represent 
them in one form 

/•OO 

T(s,t)=i bdbJ (bq) (1 - F( X (s,b))). (13) 
Jo 

In order for Eq. (I13p to produce a unitary amplitude, it is sufficient to impose 
the following properties on F: 

a) it must be an analytic function of its variable s and 6, and satisfy crossing 
symmetry, 

b) it must not exceed unity, 

c) whenx(s,6) — > 0, F — » 1, 

d) when \(s, b) — > oo, F — > 0, 

e) at large 6, it must fall faster than a power. 

The functions F(s, b) of the above unitarization schemes satisfy such conditions. 
For the eikonal it has the standard form 



F(s,b) = exp[- x (8,b)], (14) 

for the rescaled [/-matrix 

F(s,b) = 1 (15) 
1 + X{s,b) 

and for hyperbolic tangent representation 

F(s,b) = ": M - 2 f: b) hV (16) 
1 + exp(-2x(s, b) 

The differences between these functions is connected with different form of the 
saturation regime. 

In the following, we shall assume that the scattering amplitude in the impact 
parameter representation can be factorised x( s -, &) ~ h(s)f(b), and that the energy 
dependence of h(s) is a power 

h(s) ~ s A . (17) 

We find that the energy dependence of the imaginary part of the amplitude and 
hence of the total cross section depends on the form of f(b), i.e. on the s and t 
dependence of the slope of the elastic scattering amplitude. 

In Fig. 1 the process of saturation of the unitarity bound in the case of the three 
unitarization schemes is presented. It is clear that despite the slight difference in 
the growth the overlapping function with energy, the profile at the LHC will be 
similar for all these three unitarization schemes (see Fig. 2). 

The overlapping function is complex, because it must be crossing symmetric. 
In that case, s can be replaced by s exp(— i ir/2). At very low energies, before 
the first inelastic channel opens, the total cross section is equal to the elastic cross 
section and the overlapping function has a large real part. When the energy grows 
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Fig. 2. G(slhc ,b), for the eikonal (hard line and the rescaled [/-matrix (long-dashed 

line) forms at super high energies. 

and the overlapping function reaches the Black Disk Limit (BDL) (see, for example 
[3 US]); the elastic cross section will be one half of the total cross section. Then 
the real part of the overlapping function disappears at small impact parameters. 
This has a small impact on the behavior of the total cross section, as shown in Fig. 
3, where we compare the values of a tot m the case of a purely imaginary eikonal 
and of a complex eikonal. 



Fig. 3. The total cross sections of proton-proton scattering are calculated for a purely 
imaginary eikonal (plain line) and for a complex eikonal (dashed line). 



We have shown that the impact parameter is the natural variable to study uni- 
tarity of the scattering amplitude. In particular, the unitarity bound leads to the 
Black Disk Limit at high energies. The corresponding saturation phenomena may 
be connected with the saturation of the partonic density of the interacting hadrons. 
As a step in that direction, we showed that the usual unitarization schemes are in 
one-to-one correspondence with nonlinear equations. Such an approach can be used 
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to build new unitarization schemes and may also shed some light on the physical 
processes underlying the saturation regime. In the presence of a hard Pomeron the 
saturation effects can change the behavior of some features of the cross sections 
already at LHC energies. In the unitarisation schemes considered here, non-linear 
effects appear before the BDL and lead to an acceptable growth of the total cross 
sections. Saturation then leads to a relative growth of the contribution of peripheral 
interactions, and to changes in energy dependence of the differential cross sections 
for moderate values of the momentum transfer. Such effects in the differential elas- 
tic cross section can be discovered at the LHC. 
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